Theory of phonon-induced spin relaxation in laterally coupled quantum dots 
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Phonon-induced spin relaxation in coupled lateral quantum dots in the presence of spin-orbit cou- 
pling is calculated. The calculation for single dots is consistent with experiment. Spin relaxation in 
double dots at useful interdot couplings is dominated by spin hot spots that are strongly anisotropic. 
Spin hot spots are ineffective for a diagonal crystallographic orientation of the dots with a transverse 
in-plane field. This geometry is proposed for spin-based quantum information processing. 

PACS numbers: 72.25Rb, 73. 21. La, 71.70.Ej, 03.67.Lx 
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Understanding spin relaxation in coupled quantum 
dots is important for setting the efficiency of spin-based 
applications of information processing, such as spin quan- 
tum computing or controlled generation of spin en- 
tanglement 0- Phonon-induced spin relaxation has al- 
ready been studied theoreticall y in single dots for elec- 
trons 0, GL0, 0, H H GH El E3> h ° les Ei> and 

excitons Eij, and in one-dimensional coupled dots 
Recently, spin relaxation of electrons in single dots has 
been measured [l6|. 

Here we present a realistic calculation of phonon in- 
duced spin relaxation in single and coupled lateral quan- 
tum dots formed over a depleted two dimensional electron 
gas in GaAs grown along [001], the most typical growth 
direction. We show that our calculation is consistent with 
the single-dot experiment 0. We predict that: (i) Spin 
relaxation in coupled dots is strongly anisotropic with 
respect to the orientation of both an in-plane magnetic 
field (due to the interplay of the Bychkov-Rashba and 
Dresselhaus spin-orbit terms) and the dots' axis. The 
anisotropy is limited by the in-plane inversion symme- 
try only, (ii) The spin relaxation rate varies strongly 
with the inter-dot coupling, having a giant enhancement 
at a significant range of useful tunneling amplitudes due 
to spin hot spots ( anticrossings caused by spin-orbit cou- 
pling 

HE 13 This variation, which is over several 
(four to five) orders of magnitude, should be included in 
any realistic modeling of spin coherence phenomena in 
coupled quantum dots with controlled temporal evolu- 
tion of the coupling, (iii) Fortunately, the effects of (ii) 
are absent at specific configurations. The most robust 
(with respect to materials parameters) such a configura- 
tion is with dots oriented along [110] (or [110]) with the 
in-plane magnetic field along [110] ([110]). We propose 
to use this configuration for spin-based quantum infor- 
mation experiments. 

Our single-electron Hamiltonian is H = T+V + Hso + 
Hz- Here T is the operator of the kinetic energy with 
the magnetic field B introduced by minimum coupling, 
and V is the double-dot confinement potential, 

V{v) = (l/2)m^min{(r-d) 2 ,(r + d) 2 }. (1) 

The plane radius vector is r = (x,y), where x = [100] 



and y = [010] are the crystallographic axes, while d de- 
fines the distance (as well as tunneling energy at B = 0) 
and the orientation of the dots; the angle between d 
and x is denoted below as 6. The conduction electron 
mass is m and the single-dot (d = 0) confining energy is 
Hluq. The spin-orbit coupling comprises three contribu- 
tions 20]: H S o = H B r + H D + H D3 , where 

Hbr = a B R (<?xK y - a y K x ) , (2) 
H D = lc (K 2 ) (-a x K x + a y K y ) , (3) 
H D3 = (7c/2) (<J X K X K 2 - a y K y K 2 x ) + h .c ., (4) 

are the Bychkov-Rashba, linear Dresselhaus, and cubic 
Dresselhaus couplings. Kinematic wave vector operators 
are K = — iV + (e/h)A, where A is the vector potential 
to B. While both curr and the quantum average, (K 2 ), 
in the growth direction z, are tunable by a top gate, 7 C 
is a band parameter. Below we use Irr = h 2 /2maBR 
and Id = h 2 /2m"/ c {K 2 ) as effective spin-orbit lengths. 
The last term in the Hamiltonian is the Zeeman splitting 
Hz = -(i?/2)//b(t • B, expressed by the band g-factor g 
and the Bohr magneton fj,R- 

Single electron states are obtained by numerically di- 
agonalizing Hamiltonian H using the Lanczos algorithm. 
The GaAs materials parameters are used: m = 0.067m e 
(m e is the free electron mass), g = —0.44, and j c = 27.5 
eV-A 3 [23]. The linear Dresselhaus coupling is chosen to 
be j c (K 2 ) = 4.5 meV-A, corresponding to a 11 nm wide 
ground state of a triangular confining potential |2l| . The 
Bychkov-Rashba parameter a B R is 3.3 meV-A, in line 
with experiments [13,|2^]. The above ^ C {K 2 ) and cxer 
are selected to be both generic and consistent with the ex- 
periment of Ref. [13 (see Fig. EJ). O ur confining energy 
%u>o is 1.1 meV, corresponding to the confining length of 
In = 32 nm, describing the experimental system of Ref. 
[l6j. Finally, the magnetic vector potential is given in 
Landau's gauge, A = (B±/2)(— y, x, 0) for the case of 
a perpendicular magnetic field B = Bj_z; if the field is 
in plane, B = £?||(cos7, sin 7, 0), where 7 is the angle 
between the field and x , cyclotron effects are neglected. 
This is justified here for fields up to about 10 T, for which 
the magnetic length is greater than the confining length 
in the z-direction. 

While spin-orbit terms couple opposite spin states, 



2 



10" 



7-10 4 



io J 



10' 



-v — 1^ " 

1 — ■*" 

7 / 


___J — : — 1 — 1 — 1—J 1 


1 7=90° 






X Y=45 



10 15 
magnetic field B [T] 



20 



FIG. 1: Calculated spin relaxation rate for a single quantum 
dot as a function of ~B\\ applied along [110], [010], and [110]. 
The symbols are experimental data from Ref. UrA The calcu- 
lated curves for Bn > 10 T are not realistic since they do not 
incorporate cyclotron effects in the z direction. 



electron-phonon coupling enables transitions between 
such states. Here we include the most relevant terms — 
deformation and piezoelectric acoustic electron-phonon 
potentials (direct spin-dependent electron-phonon cou- 
plings appear inefficient UQ), 
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The summation is over phonon wave vectors Q = (q, Q z ) 
and polarizations A (two acoustic, t, and one longitudinal, 
I). The phonon creation and annihilation operators are 
denoted by b' and b, respectively. For GaAs the mass 
density is p — 5.3 x 10 3 kg/m 3 , the phonon velocities 
are cj = 5.3 x 10 3 m/s and c t = 2.5 x 10 3 m/s, the 
deformation potential a e = 7.0 eV, and the piezoelectric 
constant e/114 = 1.4 x 10 9 eV/m; V is the unit cell volume 
and N is the number of unit cells. The geometric factors 
M\ depend only on the direction of Q |24j . 

We calculate the rate of spin relaxation as the tran- 
sition probability (given by the Fermi golden rule) due 
to -Hep, from the upper Zeeman split ground state (de- 
noted as Tg in .19!]) to all lower states (which have nec- 
essarily opposite spin) . If the Zeeman splitting is smaller 
than the orbital excitation energy, the spin relaxes to the 
ground state (r^) only. If, however, more orbital states 
are present below the upper Zeeman split state, tran- 
sitions to all lower states contribute to spin relaxation. 
This is particularly relevant for spin relaxation in single 
dots at large magnetic fields and in coupled dots at weak 
couplings. The spin direction of a state is given by the 
sign of the expectation value of <x in the direction of B. 

In order to predict the spin relaxation rate in coupled 
dots, we first discuss a single dot case and compare it with 
experiment. This is shown in Fig. ^ where spin relax- 
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FIG. 2: (Color online) Calculated spin relaxation, in s _1 , for 
a double quantum dot along [100] as a function of a perpen- 
dicular magnetic field B± and tunneling energy /interdot dis- 
tance. The labeled contours are equirelaxational lines. Three 
spin hot spot "ridges" are visible with spin relaxation as large 
as 10 9 s _1 . The granular pattern at spin hot spots is due to 
the finite resolution of the graphics. 



ation as a function of B\ | applied at different angles 7 is 
calculated. The experiment specific value of g = —0.35 is 
taken. Spin-orbit parameters could be adjusted from a fit 
to the experimental data. However, such a fit is presently 
not possible since the calculated rates depend strongly 
on 7, reflecting the reduced symmetry (C^) of the GaAs 
interface, while the experimental data are taken for a 
single diagonal crystallographic direction, undetermined 
whether [110] or [lTO] with respect to C 2v ) These 
two directions are not equivalent, which is reflected by the 
anisotropy shown in Fig. ^ For the purposes of demon- 
stration we assume ad hoc that the experiment is done 
for [110]; the spin-orbit parameters used in this paper are 
selected (the selection is by no means unique) to quan- 
titatively describe the experiment with this 7 |26|. This 
is to demonstrate that the experiment is consistent with 
the phonon-induced spin relaxation model for reasonable 
values of spin-orbit parameters. We disagree with the 
experiment at 14 T in which cyclotron effects (beyond 
the scope of our theory) in the growth direction become 
important. 

The calculated anisotropy in Fig. appears due to 
the reduced symmetry in the presence of both Bychkov- 
Rashba and Dresselhaus terms [l^. Using a unitary 
transformation that eliminates the linear spin-orbit terms 
in a confined system 0, [2?| the Zeeman term due to B\ | 
transforms to an effective Zeeman term with magnetic 
field B e z (x,y)z along z: 



Bf = 



-B 



cos 7 sin 7 



y 



sin 7 cos 7 



(6) 



Ibr Id J \ Ibr Id 
The spin-flip probability is proportional to the square of 
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FIG. 3: (Color online) Calculated spin relaxation rate, in s" 1 , 
of a double quantum dot as a function of 7 and tunneling 
energy, for Bn — 5 T. The dots are oriented along [100]. The 
weakest relaxation is for 7 ~ 35°. Spin hot spots strongly 
influence spin relaxation at tunneling energies from 0.001 to 
0.1 meV. 
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FIG. 4: (Color online) Calculated spin relaxation, in s~\ of 
a double dot as a function of 7 and the tunneling energy, for 
B11 = 5 T. The dots are oriented along [110]. The weakest 
spin relaxation is for 7 = 135°, while at 7 = 45° spin hot 
spots do appear. 



the transition matrix element of B a z . Since the single dot 
(Fock-Darwin) upper and lower Zeeman ground states 
are coupled through the first excited orbital states which 
can be chosen to have a definite x or y symmetry, the spin 
relaxation rate is proportional to the sum of the squares 
at x and y in Eq. The spin-flip probability is then 
proportional to the inverse of the square of the effective 
anisotropic spin-orbit length £50(7), 

£sS(7) = 1/1% R + 1/& - 2 M^)/IbrId. (7) 

The period of it reflects the C^ v symmetry of the interface 
The minimum spin relaxation is at 7 = 45°, while 
the maximum is at 135°, consistent with the numerics in 
Fig. ^ The anisotropy is absent if one of the spin-orbit 
couplings dominates. Experimental observation of such 
an anisotropy would be a clear signal of a phonon-induced 
spin relaxation and could be used to extract the ratio of 
Id and Ibr- If Ibr = Id, the anisotropy is strongest — 
the spin relaxation rate due to the linear spin-orbit terms 
vanishes for 7 = 45° . Details of the analytical derivations 
will be published in a longer version of this article [2^] • 
The anisotropy of Eq. 0has been found earlier Q, while 
related anisotropics in g-factors has been predicted for 
extended two-dimensional systems |29| . 

We now move to double quantum dots described by 
the confining potential V in Eq. ^ We have already pre- 
dicted that spin hot spots in these systems appear when- 
ever the Zeeman splitting equals the tunneling energy 
(difference between symmetric and asymmetric orbital 
levels) [3. At weak coupling [d S> lo) the Zeeman split- 
ting dominates and spin relaxation proceeds through at 
least two channels, one to the symmetric (r^), the other 
to the asymmetric (To) orbital state. At large coupling 



spin relaxation is, in general, a single channel process, ex- 
cept at very large magnetic fields in which Landau levels 
form. The two regimes are separated by a spin hot spot 
in which spin relaxation is as large as orbital relaxation. 
This regime is common — typical values for the tunneling 
energy and the Zeeman splitting are of order 0.1 meV. 

Consider first double dots with a perpendicular mag- 
netic field B± which contributes both the Zeeman split- 
ting as well as cyclotron effects. The calculated spin 
relaxation as a function of tunneling energy and B± is 
shown in Fig. [5J The profile is rather complex. Spin 
relaxation is dominated by the presence of spin hot spots 
which enhance spin relaxation in most regimes of the con- 
trol parameters. For tunneling energies below 0.2 meV 
the weakest spin relaxation is for magnetic fields from 
2 to 5 T. At d = the calculated rate is that of single 
quantum dots. There is a characteristic cusp structure 
as a function of B± at a spin hot spot at B±_ «5T (such 
a cusp is not seen for the in-plane field case in Fig|T]sincc 
due to the absence of cyclotron effects the spin-hot spots 
appear at large fields of B\\ as 54 T). Our calculation 
for this single-dot case is in quantitative agreement with 
perturbative calculations 0. 

Let us now look at spin relaxation in an in-plane field 
Bu, a situation interesting for spin qubit experiments, 
since cyclotron effects are inhibited. Two important cases 
are shown for Bu = 5 T and different orientations d. The 
first case, in which the dots are along [100], is shown in 
Fig. [3J and the second case, in which the dots are along 
[110], is shown in Fig. 0] At small and large couplings, 
the spin relaxation rate is strongly anisotropic, similar to 
the single-dot case in Fig. ^ This anisotropy is greatly 
enhanced in the intermediate coupling by spin hot spots. 
In fact, spin hot spots dominate this useful regime: spin 
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relaxation is several orders of magnitude higher than in 
the single dot case (limits of either very strong or very 
weak coupling on the graph) at virtually all 7. The ex- 
ceptions are 7 ss 35° in Fig. [3]and 7 = 135° in Fig. At 
smaller (larger) Bu, the strong relaxation regime moves 
towards smaller (larger) coupling, while the two angles 
of "easy passage" remain. 

The anisotropy in both 7 and S can be explained by 
transforming the effective Zeeman field B° s , Eq. [fjl into 
the rotated coordinate system in which the x axis lies 
along d: 

Bf = -B ll x[l B 1 R cos( 1 -S)-l D 1 sm(j + S)] 

+ B\ 1 y [l^ sin( 7 - 5) - l D l cos( 7 + 5)} . (8) 

Unlike in single dots, the relevant states in double dots 
are coupled by x and y differently. At weak coupling 
the dominant term is the one containing x (which is the 
symmetry of the first excited orbital state Ta [l9|V This 
term leads to anisotropic spin hot spots and giant spin 
relaxation seen in Figs. |3and0] Spin hot spots vanish if 
' '-br cos (l — <5) — Z^ 1 sin( 7 + S) = 0, which, for 6 = 0, gives 
tan 7 = Id /I br- For our parameters the corresponding 
angle is about 35°, consistent with the numerical calcula- 
tion shown in Fig. Unfortunately, this angle depends 
on the spin-orbit coupling and is thus not robust against 
materials and growth details. On the other hand, for 
5 = 45°, the spin hot spots vanish if 7 = 135°, which is a 
universal value independent of spin-orbit coupling. This 
is confirmed numerically in Fig. 0] In the special case of 
Ibr. = Id, the condition for the weakest spin relaxation 
would be 7 = 45°, as in single dots. 

In GaAs single dots spin hot spots, which appear at 
large magnetic fields, are due to the Bychkov-Rashba 
coupling only 0,0]. On the other hand, as can be read 
from Eq. |H1 in double dots spin hot spots appear at ar- 
bitrary small magnetic fields and are caused by both the 
Bychkov-Rashba and Dresselhaus couplings [ljj, whose 
interference causes spin hot spot anisotropy. 

Similar results apply for other growth directions. In 
[111] Hd ~ Hbr |2CJ and the above results in the limit 
Id — ► 00 apply (with Ibr being a combination of both 
coupling strengths). The spin hot spots are inhibited for 
any orientation of double dots and a perpendicular in- 
plane field [cos (7 — S) = 0]; such a configuration can also 
be used in applications. In the [110] case a unique easy 
passage exists for 7 = and 5 = 90° |2s|. 

In conclusion, we have performed realistic calculations 
of phonon-induccd spin relaxation in double quantum 
dots in the presence of magnetic field. The spin relax- 
ation rate is dominated by spin hot spots in the use- 
ful regime of interdot couplings. The spin hot spot 
anisotropy allows an inhibited spin relaxation for the dots 
oriented along a diagonal of the [001] plane with a trans- 
verse in-plane magnetic field. 
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